Cellular structure in system of interacting particles. 
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The general description ol formation the cellular structure in the system of interacting particles is 
proposed. Interactions between particles are presumably well-understood and the phase transition 
in which can be studied in the scale of particle resolution. We presented analytical results of possible 
cellular structures for suspension of colloidal particles, in system particles immersed in liquid crystal 
and gravitational system. We have shown that cellular structure formation can occur in system of 
interacting particles for realistic values of temperature and particles concentration. 
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At the same time the resulting physical properties in 
the system of interacting particles have been the subject 
of much active research. The phase separation and order- 
ing of many particles system and their resulting proper- 
ties excite curiosity as fundamental physics as large num- 
ber of viable application. For description the behavior 
the many particles system with different character of in- 
teraction we must develop approach which taking into 
account a spatial nonhomogeneous distribution of parti- 
cle i.e cluster formation or cellular structure and more or- 
dering other structures. Basically, statistical description 
of many particles systems concern homogeneous states. 
The most interesting and pressing problem in condense 
matter physic is the study of phase transition with spa- 
tial nonhomogeneous distribution of particle- cluster for- 
mation. This structure was observed in usual colloids 
in the system particle which introduced in liquid 
crystal 1 and even in distribution galactic sys- 

tem Q. As evident from this examples, formation cellu- 
lar structures in the system of interacting particle have 
been large-scale demonstration for different physical situ- 
ation. Cellular structure in system particle was dramatic 
change the properties of condense matter. Depending of 
the size of particle was observed the formation of both 
reversible and irreversible clustering. During the past 
decade, colloidal system have been used as model system 
to understand phenomena the two dimensional melting 
transition. In usual colloids was observed the crystal like 
ordering, was observed hexagonal crystal structure Q- 
Experimentally, it is observed that the initial stage of 
the cellular structure originates from region with lowest 
local density of particles. In all case exist the elastic in- 
teraction leads to spatially inhomogeneous distribution 
of such particles with formation of regions of pure from 
particles. The energy interaction in system particle have 
different character and value. Minimum of free energy 
realizes cellular structure and leads to stabilization the 
regions of pure elastic media, and the particles themselves 
are placed at the boundaries of those regions. To describe 
the conditions of formation and properties of these parti- 
cle structures, one should take into account both aspects 
of the particle interaction. At first we shall give repre- 
sentation of the general approach in the description of 



behavior in system of cooperating particles. Let's show 
as in this system can be realized of a first order phase 
transition with creation of spatially non-uniform distri- 
bution of particles. The formation of different structures 
depends upon the initial concentration of particles. The 
observation of optimum cluster size has important conse- 
quences for the formation cellular structure. The system 
will try minimize the energy by breaking it apart, and 
out of many ways it can be done, the system chooses the 
nucleation and further formation of cellular structures. 

We would like show an investigating phenomena that 
occurs in colloidal system the formation of colloidal voids 
and cluster structures. Theoretical investigation of col- 
loidal system focused on thermodynamic treatment or on 
general arguments. The essential point which is neglected 
any at these arguments is that the colloidal particles are 
interacting via their direct forces and indirect interaction 
mediated by free energy of solvent phases. We propose a 
general description which include effect contains param- 
eters whose relation to formation of different structures 
in the system of interacting particles. Basically it con- 
sider very know Landau approach for first order phase 
transition coupled to the colloidal contribution and com- 
bined with the direct interaction between the colloidal 
particles. An analysis of this general model can predict a 
rich scenario of solvent induced colloidal phase separation 
with formation of cluster and cellular structures. 

The physical properties depend on the structure of on 
individual particle. However practical exist within spa- 
tial nonuniform distribution of particles, and so inter par- 
ticles interactions will be also influence the physical be- 
havior of the this distribution and also aggregation pro- 
cess by within the spatial distribution formed. Charac- 
teristic of aggregation effect very difficult, although there 
have been attempt to measure fractal dimensional form 
of distribution particles. The presentation and develop 
method for system of interacting particle, several inter- 
esting feature have become evident. The method must 
describe to observations on similar system with differ- 
ent nature of interaction between particles which formed 
systems. Based on general approach we can estimate for- 
mation cellular structure in different system. 

First from all we consider system of not interacting 
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particles. Instead of investigate the nonuniform distri- 
bution of particle we describe possibility formation the 
cavity where not particles. Empty cavity or pore in dis- 
tribution particle is very simple find and measurement. 
In the case contiguous distribution of particles with con- 
centration c the size of pore can determine as distance to 
nearest particles. We consider formation only cavity of 
spherical form. The form of cavity not play principal role 
and other form of pore we take into account after more 
precise the all process formation the cavity. For deter- 
mination of probability of formation the cavity without 
particle used the Saslaw approach [6] which can describe 
the distribution of galactic. The probability depended 
from distribution of nearest particles. For different point 
of space the probability p(r)dr of disclose the particle on 
the distance between r and r + dr is equal the proba- 
bility nonoccurrence of particles in the area with size r 
which must multiply on the probability to find the par- 
ticle on the distance r. This relation have the following 
mathematical form: 



p(r)dr 
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p(r')dr' > 4ncr 2 dr 



This equation can rewrite in equivalent form: 
-p(r) - 



d_ f p(r) 
dr \ Ancr 2 
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which have solution 
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p(r) = Aircr exp( — ) 



(1) 
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The probability p(r)dr of disclose the particle on the dis- 
tance between r and r+dr is equal the probability nonoc- 
currence of particles in volume V multiply on the prob- 
ability cdVto find the particle in the volume between V 
and V + dV: p(r)dr = cPdV. From this relation as result 
we can obtain the probability disclose the cavity without 
particle with volume V in the form 



P(V) = exp(-cV) = xp(-N) 



(4) 



for Poison distribution. The main size of area without 
particle can estimate from relation P(V) ~ h. Obvious, 
that in system of not interacting particle we can obtain 
only nonuniform distribution of particles, but we have 
non zero the probability disclose the cavity without par- 
ticle with volume V. 

Now, we can find analytical result for the probability 
disclose the cavity without particle with random selection 
the volume V in the case of system interacting particles. 
The probability the presence in system interaction parti- 
cles in ground canonical assemble N particles in volume 
V at temperature T can be write in the standard form 
0: 



W(N) = exp{(3/iN - f3F} 



(5) 



where f3 = inverse temperature, /i chemical potential, 
and F(N, V, T) is free energy which can calculate with 



the help of canonical assemble. With this probability to 
find the mean value of derived magnitude exp(qN) which 
must depend from number of particles 



(exp(qN)) =^2exp(qN)W{N) 

N 

G- 1 ex P + q)N - f3F} 

N 

= exp{ty(Pn + q) - + q)} 



(6) 



where we used the very know relation between statisti- 
cal sums of ground canonical assemble G and thermody- 
namic function: 

PV 

lnG(n,V,T) = — = *Qi,V,T) (7) 
The previous relation can rewrite in other form: 

^2z N W{N) = exp{^(zexp((3^) - V (exp(/3(j,)} = 



N 



exp{—^)exp{^ {zexp((3 ^l))) 



(8) 



Now we can expansion in series at z the right part this 
relation and to compare the equal power of z in both part 
equation we obtain 



W(N) = exp(-^)^^-(exp^ 



where 



(exp*)$ 



d(zexp(Pfi)) ' 



exp^{zexp{l3[i)) 



(9) 
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when zexp(Pfx = 0. After this operation we can obtain 
the probability find the empty cavity in distribution of 
interacting particles. This probability can write in the 
form: 



PV 

W(0) = exp(-—) 



(11) 



If we know the equation of states we can estimate in 
general case the probability existence the cavity without 
particle in the case of system interacting particles. This 
relation show that existence the cavity without particle 
depend from equation of states for particles which filled 
in this cavity. In the case nonoccurrence of interaction in 
system we have the previous result of Poison distribution. 
It is very good representative formula. If we remember 
that the Laplace pressure of bubble is P = 2jL 7 where R is 
the radius of void babble, we can obtain the probability 
in the very know form 



W(0)=exp(- 8 -^) 



(12) 



This is the probability formation the babble ow new 
phase under first order phase transition. In general case 
it is the babble of void. It is very good representation of 
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argument to correct description the formation cavity of 
void- 
First from all we consider the thermodynamic proper- 
ties of systems week interacting particle. In general case 
we can present the equation of state in the form 



PV 



b) 



(13) 



where b is the virial coefficient. The equation of state of 
the hard spheres gas can write as follows too: 



PV = kTN- 



1 



(14) 



where v 



NV 



is packing factor a particle volume is 



excluded, has quit different form though asymptotically 
is similar to above ones. In the case of system parti- 
cles which interact through short range repulsive as hard 
sphere and gravitation attractive on the far distance, in 
the limit of low packing factor, equation of state can write 
in the form : 



PV = kTN 



2 I 2Ni 



(15) 



All equations is like a Van der Waals equations of state 
with the second part of this equation representing the 
interaction in the system. If we have equations of state 
can obtain the probability existence the cavity without 
particle in the case of system interacting particles. Mean 
value of size the void or cavity without particle we can 
estimate compare probability to i. In the case nonin- 
teracting particle we can obtain cV — ln2 or number of 
particle which can put up to this void is N = ln2 < 1. In 
the system of noninteracting particles can not form the 
voids or cavity without particles. If we consider the sys- 
tem week interaction particles we can use the equation of 
states in the standard form. In this case we can estimate 
the volume of void from relation 



PV 

— =N(l-b) = ln2. 



as(^)^(l-6)=Zra2 



V = V a 



ln2 



1/(1 -&) 



(16) 
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where Vq is volume of one particles, and v as previous 
is packing factor. If v ~ 1 and b <C 1 is very small we 
can tall that in the system with week interaction can not 
formed the voids. The size of voids will be the smallest 
as size of particles. If the packing factor decrease and 
increase the energy interaction between particle to come 
into the picture the condition to formation the cavity 
without particle and the size of one cavity can evaluation 
as 



V = V Q 



ln2 



v(l-b) 



(18) 



As example, if v ~ 0, 5b or smallest and b = 0, 5 we 
obtain that V ~ 3Vb or bigger. Next example is the 
hard spheres gas. In this case of system particles which 
interact through short range repulsive as hard sphere we 
can obtain that 



V = V 



ln2(\ 



(19) 



If v — > 1 we have that V — -> 0. It is case compact packing 
hard sphere is obvious that in this system can not formed 
the voids. If v — > as result we have that V = Vb^r 
as previous result. We can assert that the formation of 
voids in the system of particles is result of interaction 
and depend from concentration of particles. In the case 
increasing of interaction of particles, when the relation 
between potential and kinetic energy is bigger and we 
can not use the description of the system in the term 
virial expansion we must used only equation of state. 

In process of aggregation can produce few cellular. We 
can estimate the number of cell in general case. If in 
process of aggregation make up m cell in the volume V 
which present the volume of pattern must hold the next 
relation: 



raV + NV = V. 
From this equation can obtain that 



V 
V 



N 



To 

V 



(20) 
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If use the previous relation between V and Vo can obtain 
the number of cell as 



l-b 2 , 
m = N——-(l 
ln2 v 



(22) 



From this equation we can conclude, that the number of 
cell without particles increase if decrease the packing fac- 
tor v and decrease the interaction energy. If packing fac- 
tor rj — > 1 that relation with crystal structures the num- 
ber of voids go to zero. If we consider the system non in- 
teracting particle bi — > we obtain that m = JVttj (1 — v) 
and depended only packing factor. If v — » 1 the number 
of cell came to zero too. The number of particles in this 
case came to zero. If 77 — >• 1 the number of cellular cam to 
zero on other motivation. In this case all particles formed 
close packing system and in this system can not formed 
the cavity without particles. 

If we have the inter particle interaction energy, we can 
study the thermodynamic behavior of an aggregate of 
such particles and describe the condition for the creation 
of new structure. The character and intensity of the inter 
particle interaction in the system of particles in different 
matter can be such that a temperature and concentration 
phase transition and produce a spatially inhomogeneous 
distribution the particles in the investigating system Q 
Q . We must describe the first-order face transition when 
the external field present surface boundary condition on 
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all particles. In order to demonstrate the mechanism 
and character of the phase transition accompanied by 
the formation of the inhomogeneous distribution of the 
system foreign spherical particles. For this case the free 
energy of particles in the self-consistent field approach 
and the many-body approximation can be written in the 
form: 

F = F p + F s + F(n) (23) 
where the first part 

F p = J U(r-?)f(f)f(f')drdr' + ... (24) 

describe the free energy which cam from interaction in 
term function f(r) spatially distribution particles. The 
simplest free energy including a interaction terms that 
respect this local gauge symmetry. Se second part 

F s = J {/(r) In f(r) + [1 - f(r)} ln[l - f(r)]dr (25) 

is the entropy part free energy. This kind the entropy 
part the free energy is motivation that the two classical 
particle not possible occupation own space place. Next, 
the free energy resulting from the coupling between dis- 
tribution function and colloidal coordinates can modelled 

by 

p n= I /(i 1 )^^--?)* (26) 
J i 

where W-^r) contains microscopical information about the 
wetting properties of the particles surface. 

The minima of the free energy corresponds to the self- 
consistent field solution for f(r). Each of thermodynamic 
functions of state corresponds to a solution that describes 
some phase of particle arrangement. If their distribu- 
tion can be inhomogeneous, then the solution serves to 
find the stable phase associated with the interaction tem- 
perature and character. The minima of the free energy 
corresponds to the self-consistent field solution for f(r). 
Each of thermodynamic functions of state corresponds 
to a solution that describes some phase of particle ar- 
rangement. If their distribution can be inhomogeneous, 
then the solution serves to find the stable phase asso- 
ciated with the interaction temperature and character. 
If the particle solution is disordered, then by definition 
the mean value f(f) = c, where c is the relative particle 
concentration. The concentration inhomogeneity gives 
rise to an additional term f(r) = c ±<£>(r) where ip (r) 
is the change of the probability distribution function of 
the particles. If the concentration inhomogeneitics arc 
smooth and their scale is much longer the inter particle 
distance, the quantity may be interpreted as the change 
of particle composition. When passing from to contin- 
uum description, we can write the free energy incre- 
ment, associated with the inhomogeneous particles distri- 
bution in the term of the power series expansion in using 



the long-wavelength expansion of the concentration i.e. 
<p (r*) = ip (r) + ~pidip (f) + \~Pi~p jdjdup (r) +. Where 

~p = ~r — r' distance between two particles. In this 
case we may be rewrite the part free energy, which is de- 
pendence from the change of the probability distribution 
function of the particles in the form: 

(27) 

Where 

,^(V-^),V = JU^W, (28) 

and 




A is the coefficient responsible for nonlincarity of system, 
which is induced the many-body interaction in system 
particles. The coefficient e — NAnR^W present the en- 
ergy which include every particle through the wetting 
effect, where Ro is size of particle and W present the an- 
choring energy of molecules of matter with surface of par- 
ticle. This coefficient can introduce if we have the same 
wetting effect on every particle. In this case in general 
presentation the free energy not exist none even terms 
because the distribution function of particles satisfy the 
relation: J f(r)df = N, J Lp(r)dr — 0. The expression 
is the Landau free energy of a system particles which 
foreign in matter, below the phase transition tempera- 
ture of this system. Thus we see that the minimum of 
the functional realizes a spatially inhomogeneous macro 
particle distribution only provided the sing satisfy some 
relation and the values of coefficients determined by the 
inter particle interaction. In order to reveal the condi- 
tion under which the homogeneous particles distribution 
become unstable, we have to calculate all the coefficient. 
Temperature of the phase transition to new states the in- 
homogeneous distribution of particles may be determined 
from following relation 

kT c = c(l - c)V (30) 

The functional (23)is very well know functional, which 
description the first -order phase transition with accom- 
panied the cluster formation in the system of the inter- 
acting particles. The most important contribution in the 
concentration is associated with the field configuration 
for which the value of the free energy is minimum, i.e.: 

dV 

A,--=0 (31) 

Where 

V = - l -ti 2 ^ + \\p i -ep (32) 
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can be used as potential energy in our case. Substitution 
the solution in the expression from the free energy yield 
and its variation due to the formation of new phase. In 
the case when the difference of minimum effective poten- 
tial values is greater than the barrier height, the solution 
in our case the free energy one cluster is described by the 
expression : 



AF = 4tt 

47T 







l(dip ) +VM 



r 3 e + 47rr 2 cr 



(33) 



Where a is the surface energy of the cluster boundary 
that is equal the free energy corresponding to the solution 
of the one-dimensional problem, i.e. 



a = / dr 




V(<p) 



tyy/W{p) (34) 



The radius of the new phase cluster by the minimum of 
the free energy. It is given by i?o = ^f- As is very know 
we can obtain e = ^^-and a = then R = ^} in our 



A3 



3A3e 



case and the effective value of the free energy variation 
from cluster formation is given by 



AF = 



8naR 2 



(35) 



They sometimes possess mixed physical properties of 
their elements, but in many cases quite new properties 
emerge, reflecting new structural organization of their 
elements. The characteristic dimensions of the spatially 
inhomogeneous distribution of the concentration of par- 
ticles in the end, become a criterion the first order phase 
transition with cluster formation in the system particle. 
The criterion of instability given by this relation can be 
interpreted as a condition for formation of spatially non- 
uniform distribution at a given temperature, which de- 
pends on the concentration of particles and characteristic 
length of the new structure. If we have the some fixed 
concentration of particles in the system the process of 
aggregation finished when all particle to assemble in few 
cluster. The size of cluster will be growth to time when 
all particle will be located only cluster. Cluster formed 
the interaction energy and distribution the particle in one 
cluster can consider homogeneous. If the concentration 
is bigger the volume of cluster growth too, and come the 
moment when the particle formed one cluster which to 
hold a position on volume of all particles. Other volume 
will be formed the voids. If the concentration is small we 
will be observer only cluster with equilibrium size, if the 
concentration is bigger as v > i we obtain the one clus- 
ter and few voids. In present approach we can estimate 
the size of voids as characteristic length of instability 
the nonuniform distribution of particle. In summary we 
have independently estimated the spontaneous formation 
loosely bound, order aggregates of colloidal particles and 
possible formation the cellular structures as result the 
different nature of interaction. 



The probability of formation one cluster can write in the 
form: 

, AF , 87rcri? 2 N , N 
P{R) = exp(- — ) = exp( — ) (36) 
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